Accurate statistics of a flexible polymer chain in shear flow. 
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We present exact and analytically accurate results for the problem of a flexible polymer chain in shear flow. 
Under such a flow the polymer tumbles, and the probability distribution of the tumbling times r of the polymer 
decays exponentially as ^ exp(— ar/ro) (where ro is the longest relaxation time). We show that for a Rouse 
chain, this nontrivial constant a can be calculated in the limit of large Weissenberg number (high shear rate) 
and is in excellent agreement with our simulation result of a ~ 0.324. We also derive exactly the distribution 
functions for the length and the orientational angles of the end-to-end vector R of the polymer. 

PACS numbers: 02.50.-r, 83.80.Rs 



Dynamics of a polymer under a shear flow has been of great 
interest both experimentally and theoretically lUl S S B B 



111 11211 . In biological systems, biomolecules 



subjected to complex fluid flows j^l [ij] are quite common, 
and a shear flow is one such example. In shear flow, a poly- 
mer gets stretched as well as tumbles in an irregular fashion. 
A crucial quantity which describes the interesting conforma- 
tional evolution of the polymer is its end-to-end vector R (see 
Fig. 1). Rec ently , experiments on a single DNA molecule 
in shear flow llOll have obtained accurate probability distri- 
bution functions of the length, the orientational angles, and 
the tumbling times of the vector R. On the other hand, the- 
oretically, although scalingresults from studies of non-linear 
single bead-spring mod el 13, 0, |^ and approximate analysis 
of semi-flexible chains ll ill are known, these are mostly non- 
exact. 

For a non-linear system as a semi-flexible polymer (like 
DNA), approximate theoretical results as in are per- 

haps as best as one can get. They agree well with the static 
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properties seen in experiments IIICJII . On the other hand, exact 
and analytically accurate results are very desirable for at least 
the flexible polymer problem. In particular, there exists no 
theory for the tumbling time statistics of the vector R; heuris- 
tic arguments given in are simply inadequate, as will 
be evident from our analysis below. In this Letter, we derive 
exact and analytically accurate results for the static and dy- 
namic properties of i2 of a flexible Rouse chain [14] in shear 
flow. 

The stochastic process of our concern, namely the end-to- 
end vector R{t) of a linear polymer is a Gaussian random 
variable and its dynamics is non-Markovian. The aspect of 
Gaussianity makes it quite easy to write down the static "joint" 
probability density function (PDF) of the Cartesian compo- 
nents of the vector R, and consequently the joint PDF of the 
length R, latitude angle 6 and the azimuthal angles |S Hj. 
The first non-triviality is to get the PDFs of the individual po- 
lar co-ordinates, namely, F{R), U{9) and Q((t>), in the sta- 
tionary state. While was known {[^[Hl], in this Letter 
we derive F{R) and U{9) exactly for a linear chain. 

Secondly, the non-Markovian evolution of R{t) ifTsIl makes 
the first-passage questions (like the tumbling time statistics) 




FIG. 1: A polymer configuration with shear flow along a::-direction. 



analytically extremely challenging. First passage questions in 
the context of polymers have been of long standing interest 
lis 12, 18]. In general, for non-Markovian processes, calcu- 
lation of first-passage properties are very nontrivial even when 
the full knowledge of the non-ex pon ential two-time correla- 
tion function is available Ellis |21L 12211 . However, when a 
process is smooth (as defined below) a method called 'Inde- 
pendent interval approximation' (IIA) is applicable 1IT9I1 and 
yields accurate estimates Il2ll22ll . Very interestingly, while in 
the absence of shear any component of R{t) is a non-smooth 
process and thus analytical prediction is unknown, we show 
below that in the presence of strong shear, a suitable com- 
ponent of R{t) associated with tumbling becomes a smooth 
process, leading to analytical tractability via IIA. Thus quite 
unexpectedly mathematical simplicity is achieved in a case of 
greater physical complexity. To be precise, we show that the 
PDF of "angular tumbling times" t, goes as ^ cxp (— ar/ro) 
(where tq is the longest relaxation time of the chain), and in 
the limit of large shear rate a — > 0.324+. This number will 
serve as a lower bound for experiments. 

As shown in Fig. 1, we study the Rouse dynamics QUI 
of a polymer chain of N beads connected by harmonic 
springs, in a shear flow in the x-direction. Let r„(t) = 



2 



[xn{t), Unit), Zn{t)]'^ denote the coordinate vector of the n* 
bead (n = 1, 2, . . . , iV) at time i. For = 2, 3, . . . , TV- 1, the 
position vectors evolve with time according to the equation of 
motion 



-T— = fc(r„+i +r„_i 
dt 



2r„) + f„{t)+r]in,t), (1) 



where k denotes the strength of the harmonic interaction 
between nearest neighbor beads, and the vector /„(^) = 
[7y„(t), 0, 0]"^ denotes the shear force field with rate 7. 
The Weissenberg number Wi = jtq, where the longest 
relaxation time tq — N'^/kn'^. The vector ri{n,t) = 
[rii{Ti,t),ri2{Ti,t),ri3{n,t)]'^ represents the thermal white 
noise with zero mean and a correlator {rii{n,t)rij{n' ,t')) = 
CSij5n.n'S{t — t'), where i,j = 1,2,3 and n,n' = 
1, 2, . . . , iV. The noise strength ( is proportional to the tem- 
perature and all the force strengths in Eq. ([T]i are scaled by vis- 
cosity. With free boundary condition, the two end-beads (for 
n = 1 and n = N) feel only one sided interaction and there- 
fore they evolve via modified equations which is obtained 
from Eq. ([TJ by using ro{t) = ri{t) and r^v+ilO — r]y{t), 
for two fictitious beads and + 1. 

For large N limit, the discrete n of the beads is replaced 
by a continuous variable s l23ll and the discrete Laplacian in 
Eq. ([T]i is replaced by a continuous second derivative along s 
direction. Equation ((TJ then leads to 



dr{s,t) d'^r{s,t) 



dt 



/(s,i) + r/(5,t), 



(2) 



with the free boundary conditions dr{s, t)/ds = at s = 
and s = N . In this continuum limit, the shear field and the 
noise correlator are given by /(s, t) = [72/(5, t), 0, 0]^ and 
{r]i{s,t)r]j{s' ,t')) = (SijS{s — s')S{t — t') respectively. The 
end-to-end vector is R{t) ~ r{N, t) — r(0, t). 

Solving Eq. dU by the Fourier cosine transformation 

oc 

v-^ , , /m'!TS\ 

r(s, t) = r(0, t)+Y, r{m, t) cos , (3) 

m— 1 

we find that the Fourier modes are given by 

f{m,t)^f f{m,t') + f}{m,t') exp[~a„,{t - t')] dt' 



(4) 



+ r(rn,0) cxp(— Omi) 



where a™ ~ tji^/tq, and f{m,t) — [7y(m, t), 0, 0]"^ and 
fi{m,t) = [i]i{in,t),7]2{m,t),r]3{'m,,t)]'^ is the m* Fourier 
mode of the noise vector ri{s, t). The zero mode f(0, t) ~ 
N^^ r(s, t) ds describes the center of mass motion of the 
polymer Equations (O and (|4]i, together with the knowledge 
of the Fourier space noise correlator {r]i{m,t)r]j{m' ,t')) = 
2C,N^^Sij6m.m'S{t — t'), after some algebra, leads to the two 
point space and time dependent correlation function between 
different relative position coordinates. In the stationary state 



limit t ^ 00 with a finite time increment r > 0, the correla- 
tion function 

- f ('■'(0,t)] • [r(j)(s2, < + T)-f(^')(0, t + r)]) 
t^oo 24 (miTSi\ nmrs2" 
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(5) 



where i,j = 1,2,3 and the notations r^^\s,t) = x{s,t), 
r^'^^s,t) = y{s,t), and r^^^{s,t) = z{s,t). Note that, any 
static correlation function in the stationary state limit can be 
obtained by setting r = in the above equation. On the other 
hand, for dynamic properties like tumbling, we need the cor- 
relators with T > 0. 

We first consider the static distribution functions re- 
lated to J? = [Rx,Ry,Rz]^- Since the Cartesian com- 
ponents are Gaussian random variables, their joint PDF 
is, r{R^,Ry,R,) = (27r)-3/2|C|-i/2cxp(-iH^C"iH), 
where C denotes the covariance matrix. Putting t = and 
(si,S2) = (0,0),(0,iV),(iV,0) and iN,N) inEq. ©we get 
the four correlators, respectively, needed to calculate each ele- 
ments of C Since the vector i? is a difference of two position 
coordinates, the zero modes cancel, and hence we need not 
worry about it being present in Eq. (|5]l. We find. 
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^^'^-wi^ (7) 



48fc 480fc 

and the determinant \ C\ = ac, with c = ad — b^. 

The scaling of {RD = w NWi^ - iV^ is a pathological 
feature of the Rouse model in shear flow in comparison to 
semi-flexible polymers (better modeled with FENE constraint 
on i2 If the Wi dependencies are ignored, the asymptotic 
functional forms of the PDF's of the Rouse model, compare 
quite well with semi-flexible polymers. 

It is straightforward to obtain the joint PDF of the polar co- 
ordinates by using the standard transformation: V{R, 9, 4>) ~ 
ViRcosO coscf), Rcosd sm4>, Rsiii9)R^ cos6. By integrat- 
ing it over R we get the joint angular PDF 
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(8) 



Now again integrating over (p in Eq. we obtain our first 
important result, the PDF of the latitude angle. 
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7r(di — d2)\/di + d2 
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di + d2 



(9) 
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where di = [a{d + a)/2] cos^ 9 + csin^ 0, d2 = ia/2)[{d - 
a)2 + 462] 1/2 cos2 e, and E(g) = /^/'(l - sin" 0) '^'d/? 
is the complete elliptic integral of second kind ll24ll . From our 
exact result in Eq. (|9]l, by taking the limit of 7 ^ 1 and 9 <^1, 
we see that the E(.) constant, and [di + ^2) ~ ad and 
(di - d2) « c9^: these lead to U{e) ~ Wi"^6'"2. Further, 
from Eq. in the same Hmit 5(</> = 0, - Wi"^6'-3. The 
azimuthal angle distribution Q{4)) = "-^^ 
easily derived from Eq. ([8]i; we skip its explicit expression as 
similar result has been derived earlier in 191 [llj]. We just note 

that Q{(j)) peaks exactly at = 0,„ = i tan^^ (^T^) ■ '^'^^ 
full width at half maximum A0 of Q{(l>) is given exactly by 
cos(2A(?!)) = 2 - [(d + a)/((i - a)] cos2(/),„ and for 7 > 1, 
A(/) « (?!)„i ~ Wi"^ and Q(^) « Wi~^ sin"^ ^he asymp- 
totic dependences of the functions U{9), S'((/) = 0, 9) and 
Q{<f>) on 9 and 0, that we derive from our exact Eqs. (O and 
^ for a Rouse polymer, match with earlier studies on semi- 
flexible polymers BSIEESIII- 

To derive our second result for the radial length distribu- 
tion function F{R) = /J'' d0 J^^^^^ d9V{R, 9, 4>), we em- 
ploy the trick of first calculating the Laplace transform of 
the PDF of i?2 instead, namely h{s) = (cxp(— si?^)) — 
(cxp(— s[i?2 -(_ _|_ /^2J~J^) jYiis is easily obtained as h{s) = 



(l + 2as) (1 + 2s(a + d) + 4cs2) The PDF 
is then related to the the inverse Laplace transform of h{s) as 
FiR) = 2i? X [£-i{/i(s)}(i?2)] and is given exactly as 



F{R) 



dx- 



Ioi^iR^x) 



where 7n(.) is zeroth order modified Bessel function of first 
ki nd B, and A = l/(2a) - (a + d)/{Ac) and = 
^{d~ aY + 462/(4c). From asymptotic analysis of Eq. (fTol i 
we see that F{R) « V2i?V V'l'IC'l fd" smafl R, and « 
exp(-[l/(2a) - A - ^]i?2)/y4^i^C|7I(^7+A) for large R. 

We now turn our attention to our third and main result on 
the first passage question of the polymer "tumbling" process 
in shear flow. The tumbling event is either defined as a radial 
return of the polymer to a coiled state (as in experiments ifioll 
and simulations j^), or as a angular return of vector i? to a 
fixed plane j^, say </> = 0. The former radial definition re- 
lies on an arbitrary choice of a threshold radius , while 
the latter angular tumbling is not. In this Letter we study 
the statistics of angular tumbling time, i.e. the distribution of 
times T between two successive zero crossings of the stochas- 
tic process Rx{i)- For the scaled time T = t/tq the relevant 
PDF asymptotically is P{T) ~ cxp(— aT). Analytical scal- 
ing dependence of a on Wi is known for semi-flexible poly- 
mers [71], but accurate constant factors were not estimated. We 
show below that for a Rouse chain, in the limit of large Wi, 
a approaches a constant value, and that can be estimated by 
using a systematic IIA calculation. 

For a Gaussian stationary process, the mean density of zero 
crossings is given by 12511 : p = — -^Z— ^"(0)/7r, where 



A(T) is the normalized correlator, i.e., ^4(0) = 1. We need 
A'{0) = and a finite A"{0) for p to be finite — then the 
process is smooth and one can use IIA ifloll . 

Now, for the relevant stochastic process R^ (t) of our con- 
cern, using Eq. (|5]l we find the stationary state correlator 
= limt^oo {Rx{t)Rx{t + roT)) as. 
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The two sums in the first and the second lines in Eq. ( fTTT ) for 
Cxx (T) will be henceforth referred to as Csh (T) (due to shear) 
and Cth (T) (due to thermal fluctuations) respectively. 

In the absence of any shear (Wi = 0), the normalized cor- 
relator becomes AiT) = C,h(T)/C,h(0). Using Cth(T) from 
Eq. (fTTT i we see that both ^'(0) and ^"(0) diverge, which in 
turn makes p infinite. Thus, in this case Rx (t) is non-smooth 
— see inset (b) of Fig. 2. Although IIA fails in this rather sim- 
ple looking case, our numerical simulation gives a ~ 1.20 
(see the curve for Wi = in Fig. 2). 

On the other hand, in shear flow (Wi 7^ 0), both the terms 
Csh{T) and Cth(r) are present in C^xiT), and the smafl T 
singular behavior of Cth(r) contribute also to Cxx{T). Thus 
although Rx{t) has long excursions (see inset (a) of Fig. 2) 
the thermal noisy contributions keep it non-smooth. While 
this makes application of IIA seem hopeless, we note that for 
(10) Wi^ ^ 2 in Eq. ( fTTT ). the term Cth can be ignored compared 
to Csh- More precisely, in the limit Wi 00, the normal- 
ized correlator A{T) Csh(r)/Csh(0). Using Csh(T) from 
Eq. one finds ^'(0) = and A"{Q) ^ -UO/n* 
a finite mean density of zero crossings p ~ \/120/ tt' 
fact that the process Rx (t) becomes smooth is clearly seen in 
the inset (c) of Fig. 2. Thus in this limit of strong shear, IIA 
becomes applicable. 

A crude estimate of a can be made by approximating 
P{T) to be exponential over the full range of T and not just 
asymptotically — this gives a ~ p = 0.353. For a more 
systematic approach one needs to use IIA. In Ref. |l9ll . 
few different IIA schemes were discussed. We calculated 
a by all these various schemes, and the various estimates 
differ slightly — these details will appear in a future publi- 
cation. In this Letter, we present a particular approximation 
which yields a very close to numerics. We start with ll26ll 
MT) ~ En=oi-^TPniT), where p„(r) is the probability 
of having n zero crossings of Rx between and T. Then IIA 
assumes pn (T) to be a product of the probabilities of intervals 
which make up the stretch to T, integrated over the locations 
of the zero crossings. The latter convolution integrals are 
best handled by Laplace transformation, and one eventually 
obtains a relation between the Laplace transforms A{s) and 
P(s), of A{T) and P{T), respectively d Q: P(s) = 
[1 - ((T)/2).(l - sA{s)]/[l + ((T)/2).(l - sA{s)]. 
From Eq. ( fTTT i. we obtain the exact Laplace transform of 
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approaches the value in Eq. (fTsT l from above. On the other 
hand, for the large Wi regime, it is known from experiments 
ifioll and FENE model simulations |0] that a increases as Wi 
increases. These two facts put together imply that a would 
reach a minimum value for some intermediate Wi and that can 
only approach the value in Eq. ( fT3] l from above. In summary, 
we have obtained exact PDFs for the length and latitude an- 
gle of the end-to-end vector of a Rouse polymer in shear flow. 
Further, we have derived an accurate analytical estimate of the 
decay constant associated with the PDF of angular tumbling 
times for a Rouse chain in the limit of strong shear. 

We thank S.N. Majumdar and A. Sain for useful discus- 
sions, and grant no. 3404 — 2 of "Indo-French Center for the 
Promotion of advanced research (IFCPAR)/CEFIPRA)". 



FIG. 2: Linear-log plot of P{T) versus T: the fitted a's for the four 
curves with Wi = 0,2.0,6.1, and 20.3 ai'e 1.20,0.57,0.375, and 
0.326 respectively. The two datasets with symbols □ and A in (c) 
are obtained by switching off the thermal noise along x direction 
(rji = 0) in Eq. (|T]( and 7 = 0.2 and 0.6 respectively — both fit well 
with the analytical a — 0.324 line. All data are for A'' = 10. Inset: 
Typical (t) versus t corresponding to the cases (a) both 7 7^ and 
771 7^ 0, (b) 7 = and i]i 7^ 0, and (c) 7 7^ and 771 = 0. 

Csh(r) and hence A{s) in the limit Wi — > 00 as. 

Since P{T) cxp(— aT), the Laplace transform P{s) must 
have a simple pole at s = — a. In other words, the denomina- 
tor of P(-a) must vanish, i.e. l-{{T)/2)a{l + aA{-a)) = 
0, where (T) = 1/p = 7r^/VT20. Solving for a from the lat- 
ter, we finally have, 

a = 0.323558.... (13) 

To check the accuracy of our analytical result Eq. ( fT3] l, we 
perform a simulation switching off the thermal noise in the 
X direction (rji = 0) in Eq. ([T]i — this effectively achieves 
the limit Wi 00 for any finite 7. For the latter case, with 
7 = 0.2 and 0.6, we show in Fig. 2 that their slopes a for 
P{T) have excellent agreement with Eq. (fTsT l. For any finite 
Wi (keeping 771 7^ 0), the value of a smoothly interpolates 
between the two limits ~ 1.20 and 0.324 (see Fig. 2). 

No direct cornparison can be made with the published ex- 
perimental data llOil , as the latter study is for radial tumbling. 
We look forward to future experiments on angular tumbling 
of a semi-flexible polymer We claim that our result for a 
in Eq. (T3[ will serve as a lower bound, based on the fol- 
lowing argument. For the small Wi regime, a semi-flexible 
polymer may be represented by the Rouse limit, for which we 
have shown (Fig. 2) that a decreases as Wi increases and 
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